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TWO SHORT METHODS FOR COMPUTING THE ORBIT 

OF A SPECTROSCOPIC BINARY STAR BY USING 

THE ALLEGHENY TABLES OF ANOMALIES. 

By F. Henroteau. 

Several methods have been proposed for deriving the elements 
of spectroscopic binary stars, but they involve either a certain 
amount of computation or a certain amount of graphical construc- 
tion, or else the use of curves that have been drawn once for all on 
tracing cloth, as in the elegant method due to King. 

Since the Allegheny Tables are used to pass from the elements 
to the computation of an ephemeris, it would seem that they 
might also be used for the reverse operation, that is, to pass from 
the radial velocities to the elements of the orbit. Two methods 
have been devised, following that idea, of which the second one 
especially is very expeditious and not only allows us to get pre- 
liminary elements but affords a criterion to test them and in a very 
short time, by successive approximations, to secure very good 
elements. The two methods are given here at length with their 
theoretical explanations. The second one has been used for 
several orbits of different eccentricities and has always given good 
results. 

First Method, for Finding the Point Corresponding 
to Periastron and the Eccentricity. 

Suppose that the velocity curve has been determined by a 
series of observations. We know that the equation of the velocity 
curve is, if the mean axis XX t is taken as axis of abscissae, 

V = K cos u, (i) 
where V is the velocity and u — v + «, v being the true anomaly, 
to the angular distance of periastron from the ascending node, K 
the semi-amplitude measured on the velocity curve with the same 
unit-length as the limes or abscissae (K is also the absolute value- of 
the ordinates of maximum or minimum, since XX S is the mean 
axis). 
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From (i) we obtain by differentiating: 

dV _. . dv . , 

_=_ jR : sl n„-;( 2 ) 

from which we have 

d/ d* A .Ksin« * w 

d7 

j— is the slope of the tangent at any point A of the velocity curve. 

It is also clear that if we first multiply -j- by — to transform v 

at ir 

into degrees, and secondly multiply dt by n oi^jj-, the mean daily 

dv 
motion (V being the period expressed in days), -=- will become the 

slope of the tangent at a point corresponding to A of the curve having 
for abscissae the mean anomalies and for ordinates the true anom- 
alies of the orbit (We shall call this the anomaly curve). Suppose P 
is the point corresponding to the periastron on the velocity curve. 
Let us call P' its projection on XX% and let A' be the projection 
of A. P'A' is the time required by the star to go from periastron 
to the position corresponding to A . So if we know P,/iX P'A ' will 
be the abscissa of the point corresponding to A in the anomaly 
curve. 

If we now take the tables for the true anomaly in elliptic orbits 
published by Dr. Schlesinger and Miss Udick in Volume 2 of the 
Allegheny Observatory Publications, for each value of the eccen- 
tricity we have the abscissae M (mean anomalies), the ordinates 
v, and the slopes of the tangents A corresponding to each value of 
M . (In other words, A, the difference between two successive true 
anomalies, differs so little from the slope of the tangent that it 
may be taken as such.) 

Let X, Xi and X 2 be the points where the velocity curve crosses 
the mean axis. If XXi < XiX if P' will be on XX X and if XX X > 
■X\X 2 , P' will be on XiX^. We may suppose P' between X and X\, 
the following reasoning being the same if it should be between X\ 
and X 2 . 

Let + { and — jj be the slopes of the tangents to the velocity 
curve at X and Xi. We have from (3) and since u at these points 
equals 270 and 90 

at X,-r= £andat Xi — — — 



d/ K l dt K 
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from which the corresponding slopes of the tangents to the anomaly 

U 
curve can be easily determined by multiplying by — as has been 

shown above. Let — Ai and A 2 be these slopes. 

If P' should be known then - n . P'X = - Mx and m . P'Xi 
= My would be the ordinates of the points corresponding to X and 
Xi on the anomaly curve (if, as is understood, the point corre- 
sponding to periastron is the origin in this curve). 

We do not know yet Mx nor My, but if we call /x . XXi = M> 
we know that Mx + My = M,. Therefore, if having computed 
Ai and A 2 we look in which column of the Allegheny Tables they 
are situated so that the sum of their corresponding mean anomalies 
Mx + My = M$ (M» being deduced easily from the velocity 
curve), we shall read at the top of the column the eccentricity, and, 
moreover, the value of Mx that we shall read from the table will 
give us the position of the periastron. It will be quite unnecessary 
for obtaining a preliminary orbit, to estimate the fractions of de- 
grees for the values of M>, Mx and My. 

The accuracy of the present method depends solely on the 
accuracy with which the tangents to the velocity curve at X and X x 
are drawn. A little difference in the way of drawing them may 
sometimes make great differences in the values of their slopes. 
A very important criterion will show us if the tangents have been 
accurately drawn, and this is that »* + »»= 180 if we call — v* 
the true anomaly at X and Vy the true anomaly at X\. We have 
indeed at X cos ( — Vx + u )= o and at X u cos (»» + «)= o and 
therefore Vy + u = 270 and — vx + w = oo°, so that by sub- 
tracting, vx + Vy = 180 . If we do not find vx + v y = 180 we 
shall have to draw the tangents more accurately until the criterion 
is satisfied. 

Example of the Computation of the Eccentricity and 
Determination of the Position of Periastron. 

From a velocity curve 

{= + 2.01 v = + 1.16 K= £■ XX t = 2 d .8o U= 5 d .o7- 

Compute the eccentricity and the position of the periastron. 
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We have 

2.01 X 25 X 5-97 =Il6 

41 X 2 X3.I416 

1.16 X 25 X 5-97 6 

2 41 X 2 X 3-1416 ' 

Mx + My = M. =^<— = 169. 
5-97 
From the Tables we obtain 

e= 0.31 
Mx= 55 vx= 89°.9 

My= 114° Vy= I4I.O 

Vx+ Vy= 23O .9 

The criterion is not satisfied, so we have to conclude that the values 
of £ and tj have not been accurately determined. The tangents at 
the curve will have to be redrawn. A little practice in the use of 
the method will soon give us very good tangents. 

Second Method For Finding the Eccentricity. 
Altho the preceding method is good it is not advisable to use it 
if the position of the periastron point can be determined without 
great error by Schwarzschild's method, for in this second method 
the use of a criterion will permit us to correct the position of that 
point and obtain very rapidly an accurate position. Schwarz- 
schild's method consists mainly in this: Lay a piece of semi-trans- 
parent paper over the velocity curve, and copy on this the curve 

and the mean axis, marking also the points o,—,U and — — . Shift- 
ing this copy bodily along the mean axis for a distance — or half the 

period, rotate the copy 180 about the mean axis — i. e., turn the 
copy face downward on the original curve, keeping the mean axis 
in coincidence, and bring the point o or U of the copy over the point 

— of the curve. The curves will then cut one another in general in 
2 

four points of which two will be the points of periastron and 
apastron. The proper two points will be found without difficulty, 
for peri- and apastron must be separated in time by one half a revolu- 
tion, and must moreover, lie on different branches of the velocity 
curve. To determine which is periastron we have the criterion, 
that the curve is for a shorter time on that side of the mean axis 
on which the point of periastron lies. 
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P'X^ = M* and P'Xx • ^ = My. 



Now, when the periastron P is determined (see figure) we have 

We then take Mx and My 

in the Allegheny Tables and find out in which column the corre- 
sponding Vz and Vy have a sum = i8o°. The eccentricity read at 
the top of the column will be the eccentricity of the orbit. More- 
over having Vy we have <o = 270 — Vy. 

Then, to finish, the equations 

y = y' - Ke cos w = y' - (e X PP') 

a sin i = [4-13833] KU V 1 - c* 
in which the constant of attraction is given by its logarithm and 
U is the period expressed in days, will give us the velocity of the 
center of mass and the product of the semi-major axis of the true 
orbit by the sine of the inclination of the orbit. 

This last method of finding the elements of a spectroscopic 
binary star is by far the simplest of any yet known. It requires a 
very small amount of computation. It is mainly a reading off the 
tables. 

For this second method we have also a very important criterion. 
If P is really the periastron point after we have w = 270 — Vy, we 

PP' 

must have cos « = -=-. If this equation is not satisfied, P is not 
A 

the periastron point, and a better position of this point must be 
found. By the use of this criterion, even without using Schwarz- 
schild's method at all, and only assuming first for P an arbitrary 
position, we may by a series of approximations find its true posi- 
tion. For instance, having computed o> = 270 — v y for the 
arbitrary position we shall have PP' = K cos w; a better position 
of P will then be that having for ordinate PP' = K cos <o; this new 
position of P will then give us a new u = 270° — Vy and again 
a new P until the true position has been found. The use of 
Schwarzschild's method to obtain a first position of P is however 
advisable. 
Lick Observatory, September, 1017. 



